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ABSTRACT 
In  this  thesis  the  generation  of  equivalent  networks  by  means  of  a 
congruent  transformation  using  a  variable  transformation  matrix  applied 
to  the  admittance  matrix  of  an  n+1  node,  common  datum,  network  is 
studied.    The  ranges  of  the  values  of  the  elements  of  the  transfor- 
mation matrix  which  yield  passive  RLC  components  for  the  equivalent 
network  without  mutual  inductances  are  considered.    When  the  trans- 
formation matrix  is  diagonal  the  solution  involves  a  set  of  linear 
algebraic  inequalities.    In  the  general  case,  when  all  the  elements 
of  the  transformation  matrix  are  present,  the  solution  involves  non- 
linear algebraic  inequalities.    The  region  of  passive  component 
realizability  consists  of  a  hyper-volume  in  the  space  of  the  elements 
of  the  transformation  matrix,  and  the  use  of  exponential  searching 
functions  to  generate  equivalent  networks  with  design  constraints 
is  investigated . 
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Equivalent  Networks  -  Historical  Background 
1 .1   Introduction 

It  is  well  known  that  there  is  no  unique  solution  to  a  synthesis  prob- 
lem for  realizing  a  given  driving  point  or  transfer  function.    A  group  of 
networks  can  always  be  constructed  which  have  the  same  input  immittance, 
or  transfer  function,  or  both.    These  groups  of  networks  are  called 
equivalent  networks  for  the  particular  network  function  of  interest.    In 
general,  equivalent  networks  can  be  generated  by  two  distinct  approaches  . 

(a)  Using  different  synthesis  procedures  for  realizing  the  network 
function,  such  as  Brune,  Darlington,  Bott-Duffin,  Miyata  or  other  pro- 
cedures for  the  case  of  input  impedance  synthesis.    These  procedures 
are  well  covered  in  the  literature  for  input  as  well  as  transfer  function 
synthesis,  Ginllemin  [3]  ,  Van  Valkenburg  [8]  . 

(b)  Using  a  proper  linear  transformation  of  the  admittance  or  im- 
pedance matrix  of  a  given  network  such  that  the  immittance  between  any 
desired  node  pair,  or  the  transfer  function  relating  the  response  in  any 
node  pair  to  the  excitation  at  another  pair,  is  kept  invariant. 

The  second  method  is  a  more  general  approach  and  is  readily  applic- 
able to  computer  design  techniques  .    It  leads  to  a  technique  for  gener- 
ating a  group  of  equivalent  networks  among  which  the  best  suited  to  the 
problem  at  hand  can  be  chosen. 

The  method  of  generating  equivalent  networks  by  linear  transfor- 
mation was  first  introduced  by  Cauer  [1]  in  1929  in  Germany  and  the 
first  English  paper  on  the  subject  was  printed  in  1931  by  Howitt  [2]  . 
Guillemin  [3]  ,  discussing  this  subject,  has  shown  that  the  same  result 
as  obtained  by  a  transformation  can  be  achieved  by  elementary  oper- 
ations on  the  rows  and  columns  of  the  network  admittance  or  impedance 
matrix.    Schoeffler  [4]  has  shown  that  the  transformation  can  also  be 
applied  to  the  branch  matrices,  and  has  also  developed  a  technique  of 
generating  continuously  equivalent  networks  [5]  . 

In  this  chapter  a  discussion  of  the  Cauer's,  Guillemin  and 
Schoeffler's  method  is  presented. 


1  .2   Equivalent  Driving  Point  and  Transfer  Impedance 

The  equilibrium  conditions  of  a  network  containing  n  +  1  nodes 

with  n  ports  defined  by  the  terminals  between  each  of  the  n  nodes  and 

datum,  can  be  written  in  the  following  matrix  form: 

[I-].-    [Y  ][V  ]  (1-1) 

where: 

[  V  ]   is  the  column  matrix  of  the  node  to  datum  voltages 

[  I  ]    is  the  column  matrix  of  source  currents  connected  between  each 

node  and  datum 
[  Y  ]   is  the  n  x  n  node  to  datum  admittance  matrix  of  network 
The  admittance  matrix  is  formed  in  the  following  manner  assuming 

no  mutual  coupling  in  the  network. 

(1)  Each  diagonal  element  y, .  is  the  sum  of  all  the  admittances 
connected  to  node  i . 

(2)  Each  off  diagonal  element  y. .    i  ^  j    is  the  negative  of  the 
admittance  connected  between  nodes  i  and  j . 

The  driving  point  impedance  of  a  network  excited  at  the  input  node 
i  only  is  given  by 

zi  =4^  (i-2) 

where: 

A       =    the  determinent  of  the  admittance  matrix  [Y  ]  . 

A..    =    the  minor  obtained  from  Aby  deleting  the  i      row  and 
the  i        column. 
Equation  (1-2)  can  be  obtained  by  expanding  (1-1)  and  using  Crammer's 
rule . 

It  has  been  shown  by  Howitt  [  2  ]   that  if  the  admittance  matrix 
[  Y  ]   undergoes  a  real  non-singular  transformation    V 

[y1]  =  [r][Y  ][r]T  d-3) 

where  [F]   is  a  square  real  non-singular  matrix  of  the  same  order  as 
[  Y  ]  and  T  denotes  transpose,  then  [  Y1]   is  the  admittance  matrix  of  a 
new  network  with  the  same  number  of  nodes  as  the  original  network. 


In  general  the  input  impedance  at  any  node  i  and  the  transfer  im- 
pedance between  any  two  ports  i  and  j  of  the  new  network  is  different  from 
that  of  the  original  network.    However  Howitt  [  2  ]  and  Schoeffler  [  4  ] 
have  shown  that  if  the  i      column  of  [  T  ]  consists  of  zeros  except  for 

the  i      term  which  is  unity^then  the  input  impedance  at  port  i  remains 

* 
invariant  under  transformation.      If  two  such  columns  are  chosen,  the 

driving  point  impedance  at  the  ports  corresponding  to  the  columns,  as 
well  as  the  transfer  impedance  between  them,  is  invariant  to  the  trans- 
formation.   Example  1-1  illustrates  the  invariance  of  the  input  impedance 
under  transformation        (1-3)  for  a  simple  3  node  circuit. 
EXAMPLE  1-1 

Consider  the  resistance  network  of  Fig.   1-la. 
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Using  the  following  transformation  matrix: 
1        a 


cr] 


0        b 


where  a  and  b  are  arbitrary  numbers  yields 


[G']   = 


2(l-a-2a   ) 


-b+2ab 


-b+2ab 


2b 


An  original  proof  based  on  matrix  partioning  is  given  in  appendix  I 


resulting  in  Fig.   1-lb  where 


b(l-2a) 
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Fig.   1-lb 


A. 


2b 


11 


in 


2b2(l-a+a2)  -  (-b+2ab)2 


(a)    Substituting  a  =  3,  b  =  -    —   will  lead  to  the  circuit  of 
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Fig.   1-lc 


(b)    Substituting  a  =  3,  b  =  -1  will  lead  to  the  circuit  of  Fig.  1-ld 
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Fig.  1  .Id 


As  can  be  seen  the  3  networks  are  equivalent,  in  that  their  input  im- 
pedance is  the  same. 

As  can  be  seen  from  example  (1-1),  the  choice  of  the  elements  of 
[  T  ]  ,  other  than  the  first  column  are  arbitrary.    If  any  of  the  elements 
in  the  transformed  network  are  zero,  the  number  of  elements  in  the  new 
network  is  different  from  the  original  network. 

If  the  network  under  consideration  contains  inductance  and  capaci- 
tances as  well  as  resistances  then  we  can  write 


[Y  ]   =    [G  ]  +    |[L_1]   +    S[C  ] 


where 

LG  ] 
[C  ] 


[L-1] 


(1-4) 


=    the  network  conductance  matrix 
=    the  network  capacitance  matrix 
=    the  network  inverse  inductance  matrix 
By  the  distribution  law  of  matrix  multiplication  [  7  ]  the  transformation 
[  T  ]   can  be  applied  to  each  of  the  above  matrices  to  obtain  the  ad- 
mittance matrix  of  the  new  network . 


LY']   =    [G']   +   |[L'_1]   if    S[C] 


Comparing  (1-5)  to  (1-4)  after  th.eJttans.form  a  ti  on  (1-3) 

[g1]  =  [r][G  ]Lr]T 
[c  ]=  cr-Kc  ][r]T 
[L''1]  =  [r][L"1][r]T 


EXAMPLE  1-2 
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Fig.  l-2a 


Consider  the  R,   L,  C  network  of  Fig.    l-2a,  with  the  input  impedance 
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Take  [T]   = 


Then: 
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;  where  2  and  -2  are  chosen  arbitrarily. 


LG-]    =  [r][G  ][r] 


[L'_1]=    [T][D  ][r]T 
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[c]     =  [r  ][c  ][r] 
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The  new  network  is  shown  in  Fig.   l-2b.  The  input  impedance  is  un- 
changed. 
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Fig.   l-2b 
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It  can  be  seen  from  example  1-2  that  the  new  network  has  one  more 
element  than  the  original.    The  circuit  configuration  and  type  of  elements 
connected  between  corresponding  nodes  is  also  different. 

Guillemin  [  3  ]  has  shown  that  an  equivalent  network  can  be  pro- 
duced by  multiplying  the  same  rows  and  columns  of  the  admittance  matrix 
by  the  same  constant.    He  has  also  shown  that  by  using  a  technique 
which  seems  to  be  similar  to  linear  programming,  an  equivalent  network 
having  all  positive  elements  may  be  produced  from  a  network  with  a 
negative  element.    The  following  example  will  illustrate  the  above 
points  . 
EXAMPLE  1-3 


Consider  the  resistive  network  of  Fig.  l-3a. 
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The  network  admittance  matrix  is 


[G  ]   = 


in  40 

If  we  multiply  the  second  row  and  second  column  by  6/5  the  new  ad- 
mittance matrix  is: 

'  60  -24  0 

[G']   =         -24  36  -12 

0  -12  20 


11 


The  new  circuit  is  shown  in  Fig.   l-3b. 
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The  following  example  will  illustrate  the  invariance  of  transfer 
impedance . 
EXAMPLE  1-4 
Consider  the  network  of  Fig.   l-4a,  connected  to  a  current  source. 
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Fig.   l-4a 


The  resistances  are  in  ohms,  capacitances  in  farads  and  inductances 
in   henries  . 
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The  transformation  matrix  [T  ]  is  arbitrarily  chosen  as: 
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Then:  ^ 


[G']  = 
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The  circuit  diagram  of  the  new  network  is  shown  in  Fig.   l-4b.    The 


transfer  impedance  has  remained  invariant. 
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Fig.   l-4b 

1 .3  Realizability 

It  has  been  shown  by  Guillemin  [  3  ]  that  if  the  transformation 
matrix  [T]  is  real  and  non-singular,  the  positive  definite  character 
of  the  stored  energy  functions  and  loss  function  of  the  network  being 
transformed  remains  unchanged."       Consequently,  under  the  condition 
that  [T]   is  real  and  non-singular,  a  passive  network  under  the  trans- 
formation of         (1-3) is  always  transformed  into  another  passive  network. 
J.D.  Schoeffler  [5   ]   has  demonstrated  that  the  transformed  network 
can  always  be  realized  by  imbedding  the  original  network  in  an  array 
of  interconnected  ideal  transformers.    However,  the  realization  of 
transformed  network  as  a  transformerless  network  with  positive  elements 
only  presents  particular  problems  and  imposes  some  restrictions  on  the 
values  of  the  elements  of  the  transformation  matrix  [  T]  .    These  re- 
strictions and  allowable  ranges  for  the  values  of  the  elements  of  [P] 
such  that  the  transformed  network  can  be  realized  with  positive  R,  L, 
C  elements  and  no  mutual  inductance  is  discussed  later. 

After  the  transformation  is  performed  and  the  admittance  matrix  of 
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the  new  network  is  obtained,  it  is  possible  to  determine  whether  the 
new  network  can  be  realized  by  positive  R,  L,  C  elements.    The  follow- 
ing are  the  necessary  and  sufficient  conditions  for  [  Y' ]  to  be  realiz- 
able with  positive  R,  L,  C  elements  ,  Howitt  [  2  ]  . 

(1)    The  off  diagonal  elements  of  each  of  the  matrices  [G1]  , 
[  C  ]  and  [  L'      ]   has  to  be  less  than  or  equal  to  zero,  i.e . : 


ij 


^     0  for  all  i  ^  j 


G'..        ^0  for  all  i  ?  j 

U 

L'"1  <     0  for  all  i  ^  j 


(l-7a) 

(l-7b) 
(l-7c) 


(2)    The  sum  of  all  the  rows  of  each  of  the  matrices  [G  ]   , 


[  C  ]  and  [  L      ]   has  to  be  non-negative ,  i.e.: 


n 


z.     g'. 
j=i 


ij 


(l-8a) 
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©' 


1] 


(l-8b) 
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ij 
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U-8c) 


Conditions  (1-7)  and  (1-8)  can  be  checked  by  a  computer  after  the 

transformation  is  performed. 

1  .4   Transformation  of  the  Branch  Matrices 

Since  the  elements  of  a  network  are  explicitly  presented  in  the 
branch  admittance  or  impedance  matrices,  it  is  sometimes  advantageous 
to  perform  the  transformation  on  the  branch  matrices  . 

The  branch  current  and  voltage  matrices  of  a  network  can  be 
written  as: 


I 


bl 


V 


bl 


\m 


!bi 


V,    = 


V 


bi 


bn 


V 


bn 
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where  .1. and  V  .represent  the  current  through  and  the  voltage  across 
branch  i  respectively.    The  equilibrium  conditions  of  the  network  can 
then  be  written  as 

[Vb]  =    [Zb][lb]  (l-6a) 


or 


[Ib]    =    [Yb][Vb] 


(l-6b) 


where  [  Z.  ]   is  the  branch  impedance  matrix  formed  as  follows: 
b 

Each  diagonal  element  Z..  is  the  total  impedance  in  branch  i. 
Each  off  diagonal  element  Z..  is  the  mutual  inductance  common  to 
branches  i  and  j  and 


b  b 


-1 


(1-7) 


If  the  loop  current  matrix  of  the  network  is  defined  by 

I. 


[I  ]  = 


m 

where  I.  is  the  current  in  loop  i  then  [I  ]  and  [  I,  ]  are  related  by 
l  b 


Clb]  =    CarCl  ] 

where  [  tx  '    is  the  network  tie  set  matrix  formed  as  follows. 

+1 


(1-8) 


Each  element  a 


U 


if  I.  .  and  I.  are  in  the  same  direction 
bi  i 

if  I.  ,  and  I.  are  in  the  opposite  direction 
bi  i 


0         if  branch  i  is  not  included  in  loop  i 
Similarly  if  the  cut  set  voltage  matrix  [  V  ]   of  a  network  is  de- 
fined by 


V 


1 


[V  ]    = 


V. 


V, 
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where  V.  is  a  voltage  with  an  arbitrary  direction  assigned  to  the  branch  i 

of  a  tree  of  a  network,  then 

[vb  3=  L£]T  [V  ]  (1-9) 

where  [  jS  ]   is  the  network  cut  set  matrix  formed  as  follows: 

+  1  if  V,  .  and  V.  are  in  the  same  direction 
bi  i 

Each  element  /3. .    =^-1  if  V,  .  and  V    are  in  the  opposite  direction 

bi  i 

0  if  branch  i  is  another  branch  of  the  tree 

which  includes  V. 
i 

Schoeffler  [  4  ]has  shown  that  if  a  matrix  [  E  ]   is  formed  such  that 

cn]  =  [p]T  [r][a]  d-io) 

where  [I  ]   is  the  transformation  matrix  defined  in  section  1.2,  then 

instead  of  applying  [  T  ]   to  the  network  admittance  matrix,  we  can 

apply  [  n  ]   to  the  network  branch  admittance  matrix  [  Y,  ]  . 

b 

The  branch  admittance  matrix  [Y'    ]   of  the  equivalent  network  is 
then  given  by 

[Y'b]  =  [n][Yb][n]T  d-ii) 

The  advantage  of  this  method  is  that  the  elements  of  the  net- 
works are  readily  seen  from  [Y1.  ]  and  the  need  for  transformers  in 

b 

the  synthesis  of  the  equivalent  network  is  apparent  from  presence  of 
off  diagonal  elements  in  [Z'    ]  .    The  disadvantage  lies  in  the  fact  that 
[  n]   has  to  be  formed  from  equation  (1-10)  which  means  that  cut  set  and 
tie  set  matrices  have  to  be  constructed  before  [E  ]   can  be  formed. 
1  .5    Continuously  Equivalent  Networks 

J.D.   Schoeffler  [5  ]   has  developed  a  method  whereby  a  set  of 
equivalent  networks  is  produced  as  a  function  of  an  independent 
variable  x.    The  elements  of  each  member  of  the  set  of  equivalent  net- 
works are  functions  of  the  independent  variable  x.    For  any  value  of  this 
variable  a  network  equivalent  to  the  original  network  is  obtained. 

Consider  that  the  original  [Y  ]  matrix  is  Y     .    Define  a  new  matrix 

o 

Y(x)  whose  elements  are  some  functions  of  x  such  that 

Y(0)    =    Y 
o 
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Choose  Tsuch  that:    (where  [  U  ]  is  the  unit  matrix) 

[F]   =    [U  ]+  [B  ]Ax  (1-12) 

now 

[Y'(x)]   =    [F][Y(x)][r]T 

=  ([U  ]+  [B  ]Ax)  [Y(x)]([U  ]+  [B  ]Ax)  T 

=  [Y(x)]+  ([Y(x)][B  ]T  +  [B  ][Y(x)])Ax 

+  LB  ][Y(x)   ][B  ]TAx2 

If.  Ax    is  small  then[Y'(x)]will  differ,from[Y(x)]by  a  small  amount. 
Define: 

lAY(x)]  =  [Y'(x)]-[Y(x)]  (1-13) 

then  Equation  (1-10)  becomes: 

LAY(x)] 
Ax        =    [B  ][Y(x)]+  [Y(x)][B  ]  +  [B  ][Y(x)][B  ]     AxZ 

as    Ax  — ►      0 


to(x)] 

-£ =    LB  ][Y(x)]+  [Y(x)][B  ]'  (1-14) 

If  the  elements  of  B  are  chosen  arbitrarily  Equation  (1-11)  can  be 
integrated  to  yield  the  elements  of[Y(x)]as  a  function  of  x.     Schoeffler 
integrates  Equation  (1-11)  numerically  and  changes  the  assumed  values 
of  the  elements  of  B  as  the  integration  proceeds.     In  the  numerical  in- 
tegration process  the  values  of  the  elements  of  B  are  chosen  so  that 
the  new  matrix  Y'(x)  is  dominant.*    Thus  the  elements  of  B  are  actually 
a  point  function  of  x  as  the  integration  process  proceeds. 
For  a  network  composed  of  R,  L  and  C  we  can  write 

[Y(x)]  =    LG(x)]+S[C(x)]+|[L'1(x)]  (1-15) 

* 
i.e.  the  realizability  conditions  of  Sec.   1.3  are  satisfied. 
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substituting  Equation  (1-12)  in  (1-11)  then  expanding  and  equating 
terms  of  equal  powers  of  S  on  both  sides  we  obtain: 

d[^x)]    =    [B  ][G(x)  ]   +    [G(x)  ][B  ]T  (l-16a) 

d[^x(x)]    =    [B  ][C(x)   ]   +    [C(x)  ][B  ]T  (1-1 6b) 

[L^X(X)]=    [B][L~aiix)  ]  +    [L_1(x)][B]T  (l-16c) 


which  can  be  shown  to  reduce  to  the  form 

d[^xX)  ]  =    CM  ][G(x)  ]  (l-17a) 

d[^xX)  ]  =    LM  ][C(x)  ]  (l-17b) 

d[^X)]=    [M3CL-V)  ]  (1-170 


or 


d[^xX)   ]     =    LM   ][Y(x)  ]  (1-18) 


where  [  M  ]   is  a  matrix  whose  elements  are  some  combination  of  the 
elements  of  [  B  ]  .     Matrix  [  M  ]   can  be  found  by  solving  the  matrix 
equations: 

[B  ][G(x)  ]   +    [G(x)   ][B  ]T    =    [M  ][G(x)   ]  (1-19) 


or 


i.e. 


[B][Y(x)]+    [Y(x)][B]T      =    [M][Y(x)]  (1-20) 


[M  ]   =    [B  ]   +    [Y(x)   ][B  ]T  [Y_1(x)   ]  (1-21) 

Substituting  Y(0)  =  Y    in  Equation  (1-21)  yields  matrix  [M  ]  which 
can  then  be  substituted  in  Equation  (1-18).    After  integration  this 
yields  the  admittance  matrix  [  Y(x)  ]   of  the  equivalent  networks  as 
function  of  the  independent  variable  x. 
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The  Diagonal  Variable  Transformation  and  Equivalent 
Networks  with  Constraints 
2  . 1   Introduction 

In  generating  equivalent  networks  using  Equation(  l-3^)or  in  gen- 
erating continuously  equivalent  networks  by  integrating  Equation (l - 1 5, ) 
the  designer  has  no  criterion  for  choosing  the  elements  of  the  transfor- 
mation matrix  [  P  ]  to  achieve  a  design  specification  such  as  changing 
the  value  of  an  element  connected  between  two  nodes  in  a  given  manner 
or  eliminating  a  given  element  from  the  network.    The  designer  is  not 
even  able  to  choose  the  elements  of  [r]   such  that  the  transformed 
network  can  be  realized  by  positive  R,  L,  C  elements  with  no  mutual  in- 
ductance.   In  this  chapter  a  simplified  transformation  matrix,  namely  a 
diagonal  matrix,  is  chosen  for  the  transformation.    The  acceptable  range 
of  elements  of  [r]   such  that  the  transformed  network  can  be  realized 

by  positive  R,  L,  C  elements  only  is  found.    By  choosing  the  elements 

A  x 
of  [r]  to  be  exponential  functions  of  a  real  variable  x,  i.e.  e        ,  a 

group  of  equivalent  networks  are  generated  and  it  is  shown  that  by 

choosing  the  values  of  A  's  certain  constraints  can  be  imposed  on  the 

network.    However  as  we  shall  see  in  Chapter  4,  the  diagonal  form  can't 

always  generate  equivalent  networks. 

2  .2   The  Diagonal  Transfonra  tion  and  Realizability  Conditions 

The  transformation  matrix  is  considered  to  be  a  diagonal  matrix  with 

positive  elements,  i.e. 


11 


0 


[r]=  [rr  = 


r22 


(2-1) 


nn 
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where 


V..       >        0  for    i    =    1   .    .    .    .  n  (2-2) 


The  admittance  matrix  of  the  transformed  network  then  becomes: 

[Y1]  =    [T3LY  ][T]  (2-3) 

In  general 

[Y  ]   =    [G  ]+S[C  ]  +  |[L_1]  (2-4) 

Substituting  (2-4)  in  (2^3)  we  obtain 

[Y']=[r][G][r]+s[r][c][r]  +  |[r][L"1][r] 

(2-5) 
Since  in  general  [  Y1  ]  will  also  be  composed  of  R,  L  and  C  then 
we  can  write 

[Y']   =    [G'  ]+S[C']+|[L'_1]  (2-6) 

Comparing  (2-6)  and  (2-5)  we  can  write: 

LG']   =    [T][G  ][T]  (2-7a) 

[C1]   =    [I?][,C  Ur]  (2-7b) 

[L'_1]=  [r][L_1][r]  (2-7c) 
A  typical  term  of  transformed  matrices  [G1]  ,  [C1]  and  [  L'      ]  is 
of  the  form 

g'. .    =  y..  y. .  g..  (2-8) 

By  assumption  of  Equation(2-2)all  the  y  's  are  positive.     Hence 
the  sign  of  the  elements  of  the  transformed  matrices  [G1  ]  ,   [  C']and 
[  L'      ]  remain  unchanged.    If  the  original  network  was  made  up  of 
passive  elem^pts  the  off  diagonal  elements  of  the  [  G  ]  ,   [  C  ]  and 
[L      ]   matrices  were  negative.    Hence  the  off  diagonal  elements  of 
the  transformed  matrices  are  also  negative  and  one  of  the  realizability 
conditions  stated  in  section  1.4  namely  Equation (l-6)is  always  sat- 
isfied and  the  circuit  elements  connected  between  any  two  nodes  other 
than  node  to  datum,  remain     positive  after  the  transformation. 

The  sum  of  the  i      row  of  the  transformed  matrices  [  G'  ]  ,   [  C'  ] 
and  [L'      ]  are  given  by  the  following  equations: 
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n  n 

S'.(G)     =     S  g\.     =     E    y.  y. .  g..  (2-9a) 

j=1    13  j=1    11    jj     U 

n  n 

S'.(C)     =    EC'..        =    2     y.  y..  C.         (2-9b) 
j=1     u  j=1      n    33      13 

-1         n    -1  n  -1 

S'.(L     )=    £L..  =    E     y.  y.  L. .  Ifi    Q  . 

i  j=1  13  .    1      ii    JJ     13  (2-9c) 

It  is  clear  from  the  method  of  the  construction  of  the  admittance 

matrix  that  the  above  sum  of  the  rows  corresponds  to  the  value  of  the 

element  connected  between  node  i  and  datum,  i.e.  g'.    ,  C.     and 
,  10         10 

L'  .       . 
10 

If  the  second  of  the  realizability  conditions  that  is  Equation  (1-7) 

is  to  be  satisfied  we  must  have: 

n 

S'.(G)     =    I    y.y.g..        ^     0  for  all  i        (2-1  Oa) 

i  'n  '))  "13 


n 

S'.(C)     =    I    y.y.C        >       0  for  all  i         (2-10b) 

1  '11  '))      13 


-1  n  -1 

S'.(L     )    =     L    y.y.  .L.       ^     0  for  all  i        (2 -10c) 

1  11  '))      13 

Equations (2- 10)  represent  a  set  of  constraints  on  the  elements  of  the 
transformation  matrix.     Since;from  Equation  (2-2)  y  is  always  positive 
it  can  be  cancelled  from  the  set  of  Equations  (2-10).     Thus  in  general 
there  are  3n  linear  equations  that  have  to  be  satisfied,  n  equations  for 
each  of  the  [G1   ],   [C   ]  and  [  L'      ]  matrices  comprising  [Y1].    If 
Equations  (2-10)  are  satisfied  with  the  equal  sign  then  the  corresponding 
element  connected  between  the  corresponding  node  and  datum  is  zero. 
For  example,  if  S'    (G)    =    g    _    is  zero  then  the  conductance  connected 
between  node    1  and  datum  is  eliminated  from  the  circuit.     From  the 
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above  we  can  conclude  that  there  are  no  set  of  y  's  which  would  satisfy 
all  of  the  3,n  equations  of  (2-10)  with  the  equality  sign,  since  this 
will  result  in  an  equivalent  network  with  no  node  to  datum  elements.  , 
i.e.  with  infinite  input  admittance  at  any  port. 

The  3n  inequalities  of  the  form  (2-10)  together  with  the  n  in- 
equalities of  the  form  (2-2)  are  a  set  of  4n  inequalities  which  must  be 
satisfied  by  any  set  of  y's  to  assure  the  fact  that  the  transformed  net- 
work is  realizable  with  passive  elements  only.     Obviously  these  in- 
equalities are  all  satisfied  if  all  of  the  y  's  are  unity,  since  the  original 
networks  results  in  this  case. 

From  the  discussion  on  Chapter  1  it  is  obvious  that  to  keep  the  in- 
put imittance  between  node  k  and  datum  constant  y    must  be  equal  to 

K. 

unity  and  to  keep  the  transfer  impedance  between  parts  k  andi.  con- 
stant we  must  choose  y    and  y    to  be  unity.    There  are  then  n-1  or  n-2 
remaining  y  's  to  be  chosen  satisfying  Equations  (2-2)  and  (2-10).    For 
example,  in  the  case  of  transfer  function  equivalence  there  are  a  total 
of  (n-2)  +  3n  =  4n  -  2  inequalities  to  be  satisfied  by  n-2  unknown 
y  's.    These  inequalities  are  considered  in  the  next  section  to  determine 
if  a  solution  is  possible,  and  if  so  to  ascertain  specific  numerical  ranges 
for  the  y  's  over  which  the  inequalities  are  satisfied.    These  ranges 
define  the  domain  of  the  set  oi  y's  for  which  the  transformation  of  (2-1) 
generates  networks  with  passive  components  for  which  the  transfer 
impedance  between  ports  i  and  j  remains  invariant.    Obviously  a 
solution  cannot  always  be  expected  since  the  transformed  network  is 
confined  to  the  class  obtainable  by  the  diagonal  transformation  matrix 
(2-1).    The  general  case  where  a  complete  n  x  n  transformation  matrix 
is  used  will  be  discussed  in  Chapter  3. 
2  .3   The  Domain  of  the  Solutions 

Each  of  3n  inequalities  of  Equation  (2-10)  represents  a  linear 
constraint  on  (n-1)  or  (n-2)  y  's  depending  on  whether  the  input  im- 
pedance or  transfer  impedance  is  to  remain  invariant  under  the  trans- 
formation.   When  any  of  the  inequalities  (2-10)  are  satisfied  with  the 
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equality  sign  the  resulting  equation  represents  a  hyperplane  in  (n-2) 
dimensional  y  -  space.*    This  hyperplane  divides  the  (n-2)  dimensional 
y  -  space  into  two  distinct  regions  .    The  points  on  one  side  of  the 
hyperplane  satisfy  the  corresponding  inequality.    For  example,  if  the 
inequality  g'.      >  0    is  used  the  equation  of  the  hyperplane  is  g'.      =0. 

LQ  LQ 

Points  on  one  side  of  the  plane  yield  g'.     greater  than  zero  and  points  on 

the  other  side  yield  g'       less  than  zero.When  all  of  the  inequalities 

10 

(2-10)  are  considered ,     the  domain  of  the  allowable  set  of  y  's  is  de- 
termined by  a  hypervolume  in  y-space  bounded  by  a  set  of  at  most 
(4n-2)  hyperplanes  ,  3n  of  which  are  of  the  form; 

n 

I     y.  g..    =    0  (2-11) 

j=l      n     1J 

and  the  remaining  number  in  the  form 

y.    =    0  (2-12) 

'11 

where  i  f  k  and  i  ^  I   since  y       -  y      =1    due  to  the  required  in- 

variance  of  the  transfer  function.    Points  on  the  planes  representing 

Equation (2 -12)  are  not  acceptable  since  these  points  violate  the  required 

non-singularity  of  the  transformation  matrix  (2-1). 

Any  point  on  the  hypervolume  described  above  represents  a  set  of 
y  's  which  would  satisfy  all  of  the  inequalities  of  (2-2)  and  (2-10). 

The  realizability  volume  may  be  determined  in  a  systematic  manner 
as  follows  .     Solve  all  the  possible  combinations  of  the  (n-2)  simultaneous 
equations,  taken  from  the  set  of  (4n-2)  equations  of  the  form  (2-11)  and 
(2-12)  for  a  set  of  (n-2)  y  's,  and  check  whether  these  values  also 
satisfy  the  remaining  (4n-2)  -  (n-2)  =    3n  inequalities  which  are  not 
used  in  choosing  the  (n-2)  simultaneous  equations.    Any  set  of  y  's 
which  satisfy  these  conditions  is  a  vertex  of  the  realizability  volume. 
The  hypervolume  is  bounded  on  the  sides  by  any  of  the  hyper  surfaces 

* 

For  convenience  of  writing,  the  case  of  transfer  impedance  equi- 
valence is  considered.     The  case  of  input  impedance  equivalence  follows 
in  the  same  manner. 
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(2-13) 


which  passes  through  (n-2)  vertices. 
In  general  a  total  of 

(4n  -  2)  ! (4n  -  2)! 

[  (4n  -  2)  -  (n  -  2)]  !    (n  -  2)   !         3n!  (n  -  2) ! 

linear  simultaneous  equations  must  be  solved  and  checked  in  the  above 
manner.    In  practice  many  of  the  solutions  are  trivial  and  the  process 
may  be  programmed  on  a  digital  computer.    After  all  vertices  of  the 
realizability  hyper-volume  are  determined,  the  boundary  hyper-planes  may 
be  determined  from  the  simultaneous  equations  whose  intersections 
determined  the  vertices.    The  point  where  all  of  the  v  's  are  equal  to 
one  always  satisfies  the  realizability  conditions  and  is  always  in- 
cluded in  the  realizability  hyper-volume. 

If  there  is  no  set  of  y  's  which  would  satisfy  the  Equations  (2-2) 
and  (2-10)  we  can  deduce  that  the  transformation  of  (2-1)  will  not  pro- 
duce any  equivalent  network  with  passive  components,  and  the  real- 
izability hyper-volume  in  this  case  is  reduced  to  a  single  point  where 
all  of  the  y  's  are  unity.    If  a  realizability  volume  does  exist  it  is  possible 
to  remove  at  least  one  node  to  datum  element  by  choosing  a  set  of  y  's 
represented  by  a  point  on  one  of  the  boundary  hyper-planes  of  the  form 
(2-11). 
EXAMPLE  2.1 


Consider  the  RLC  network  shown  in  Fig.  2.1, 


M. 


m 


i 


i/2&  i 


i 


(4) 


Fig.  2.1 

where  all  the  conductances  are  in  mhos,  all  the  capacitances  in  farads 
and  all  the  inductances  in  henries. 
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[G  ]   = 


2 

-1 

0 

0 

1 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

[c  ]  = 


2 

0 

0 

0 

0 

1 

-1 

0 

0 

-1 

3 

:o 

0 

0 

0 

2 

[L"1]  = 


0 

-1 

4 
-1 


0 
0 

-1 

1 


It  is  desired  to  keep  the  transfer  impedance  invariant  so  let: 


[r]  = 


The  transformed  matrices  [G1],   [C']and  [  L'      ]are: 


1 

0 

0 

0 

0 

y22 

0 

0 

0 

0 

y33 

0 

0 

0 

0 

1 

[G']  = 


2 

~y22 

0 

0 

y22 

y       2 

722 

0 

0 

0 

0 

2 
V33 

0 

0 

0 

0 

1 

[C1   ]  = 


y22 

V22y33 

0 

~V22y33 

3^332 

0 
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[L"X]  = 


y 
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"y22V33 


">'22y33         ° 


4y 


33 


-y 


33 


-y 


33 


The  inequalities  (2-10)   in  this  case  are  12  inequalities  but  only  8  of 
them  contain  the  variables  y    and  y    and  of  these  only  6  are  distinct  and 

Li  O 

are  tabulated  below. 

From  [  G'  ]   we  have 


(a)  2  -  y^       *    0 

(b)  -1  +y22    £     0 
From  [  C '  ]  we  have 

^    0 

k  o 


(C)    V22-y33 

(d)    -y22  +  3>/33 
-1 


0 


From  [  L'      ]  we  have 

(e)  -y22+4y33"1       * 

(f)  -v33    +  1     >    0 

In  the  v-plane  the  above  inequalities  define  regions  as  shown  in  Fig.  2.2. 
Taking  into  account  the  Equation(2-2)which  states 


y22    >0 


y33    >0 


the  realizability  hypervolume  in  this  two  dimensional  example  will  be  a 
four  sided  polygon  shown  shaded  in  Fig.  2.2. 

To  solve  the  above  problem  with  a  digital  computer  we  must  solve 
all  the  combinations  of  the  six  equations  obtained  from  the  inequalities 
a ,  b ,  c ,  d ,  e  and  f  two  at  a  time . 

Table  2.1  presents  all  the  possible  solutions  of  the  group  of  two 
simultaneous  equations  together  with  the  letter  representing  the  solutions 
on  the  y-plane  of  Fig.  2.2.    The  solutions  which  yield  valid  inter- 
sections which  satisfy  all  of  the  other  inequalities  are  indicated  by  the 
asterisks  and  are  the  vertices  of  the  realizability  area.    The  boundaries 
of  the  volume  are  given  by  the  lines  representing  a,  b,  e  and  f. 
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TABLE  2 . 1 


EQUATIONS 

PAIR 

y22 

a  b 

none 

a  c 

2 

a  d 

2 

a  e* 

2 

a  f* 

2 

b  c* 

1 

b  d 

1 

b  e* 

1 

b  f* 

1 

c  d 

0 

c  e 

1/3 

c  f* 

1 

d  e 

3 

d  f 

3 

e  f 

3 

SOLUTIONS 


33 


none 


2 

2/3 

3/4 


/3 

/2 


/3 


INTERSECTION 
POINT 


none 
A 
B 

C* 
D* 
E* 
F 

G* 
E* 
0 
H 
E* 
K 
K 
K 


27 


*r 
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o 

Q. 

i 

0> 


O 
N 
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2.4  The  Exponential  Searching  Function 

If  a  volume  of  realizability  does  exist  then  for  each  point  in  this 
volume  there  exists  an  equivalent  network.    These  equivalent  networks 
can  conveniently  be  investigated  by  letting  the  elements  of  our  trans- 
formation matrix  [  T  ]  to  be  exponential  functions  of  an  independent 
real  variable    x.    Thus  let 


y. . 
11 


X,.x 

=    en 


(2-14) 


where  X..    is  a  constant  which  may  be  chosen  arbitrary  or  in  accordance 
with  some  additional  constraint  placed  upon  the  desired  network. 

The  exponential  function  of  Equation  (2-14)  satisfies  the  require- 
ment of  Equation  (2-2)  in  that  it  is  greater  than  zero  for  all  finite 
values  of  the  independent  real  variable  x.    Hence  the  transformation 
matrix  [T]  remains  real  and  non-singular  for  all  finite  values  of  x. 

Moreover  when  x  =  0,  y .  =  1 ,  and  for  this  value  of  the  variable  x  the 

n 

original  network  is  reproduced. 

When  the  exponential  functions  of  Equation  (2-14)  are  used  as 
elements  of  the  transformation  matrix,  [T]  becomes 


[r(x)]  = 


X.  ,x 
e   11 


0 
3L2 


\nn 


(2-15) 


,Annx 


A  typical  element  connected  between  two  nodes  in  the  transformed 
network  becomes: 


ij 


(X.. 

+ 

X..)x 

li 

jj 

q. . 

e 

u 

(2-16) 
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And  a  typical  node  to  datum  element  becomes: 


n 


n 


g1 .      =    I     g, .    =    L    g. .    e 
10      i-0     1J       j=l    1J 


(  X..  +X..)x 
11        jj 


(2-17) 


And  a  typical  equation  of  the  hyperplane  used  in  determining  the  real- 
izability  volume  becomes: 


n 

i 

j=i 


gue 


X..  x 
11 


=    0 


(2-18) 


g'       and  g1.      can  easily  be  plotted  by  a  digital  computer  as  a  function 
ij  10 

of  x.    The  range  of  x  for  which  all  of  the  3n  equations  of  the  form  of 
Equation  (2-17)  remain  non  zero  is  the  range  for  which  transformation 
of  (2-15)  yields  equivalent  networks  with  positive  R,  L,  C  elements. 
Differentiating  (2-15)  with  respect  to  x  yields: 


d  [T(x)] 
dx 


=   cr][x]  =  [xnn 


(2-19) 


where: 


[X]  = 


11 


L 


0 


22 


0 


nn 


Differentiating  (2-3)  with  respect  to  x  yields 

d[Y'] 


dx 


-  LY  ][!•]+  [HLY]-*1 


dx 


dx 


(2-20) 


(2-21) 
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Substituting  (2-19)  in  (2-21)  we  obtain: 
d[Y'] 


dx 


=    LX][Y']   +    [Y'][X]  (2-22) 


Thus  the  use  of  the  exponential  search  function  for  generating  equivalent 
networks  is  the  same  as  integrating  Equation  (2-22)  with  the  initial 
condition  that  [  Y'  ]  =  [  Y  ]  when  x  =  0  . 

Comparing  Equation  (2-22)  with  Equation  (1-11)  we  observe  that 
the  two  are  identical  with  [  B  ]  replaced  by  [  X  ]  . 

The  choice  of  the  X  's  may  be  made  to  impose  additional  constraints 
upon  the  transformed  network  as  follows. 

(a)  To  keep  the  transfer  impedance  between  port  k  and  the  output 
voltage  between  nodes  i  and  m  constant  (no  common  node 

between  input  and  output)  set  X,  ,    -X-.  ~   X        =0. 

kk         11  nn 

(b)  To  keep  the  admittance  between  nodes  i  and  j  invariant  under 

the  transformation  set  X..  =  -X..   .    From  Equation  (2-16)  it 

ii  JJ 

can  be  seen  that  in  this  case       g', .  =  g. . . 

U         ij 

(c)  From  Equation  (2-16) 

(X..  +  X..)  x 

g\.    =   g..  e       n        JJ  (2-23) 

ij  U 


and 


«■«  "  \ieUkk  +  Xu)X  (2-24» 


Dividing  (2-23)  to  (2-24)  we  obtain 

[(x..+x..)-(xkk  +  x^] 


X 


9'ij/9'k£  =  9ij/gk  e 

Hence  to  keep  the  ratio  of  the  admittance  connected  between 

nodes  i  and  j  to  the  admittance  connected  between  nodes  k 

and  &  constant  set  (  X. .  +  X. .)  =  (  X,  ,    +  X„J  . 

ii         jj  kk         JUL 
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(d)    Since  the  circuit  elements  between  nodes  i  and  j  increase  or 

decrease  monotonically  with  x  depending  upon  the  sum 

(A..  +  X..  ),  the  X  's  can  be  chosen  so  that  a  specified  com- 
11        JJ 

ponent  changes  more  rapidly  or  slowly  compared  with  others 

over  the  range  of  x. 
The  use  of  the  exponential  search  function  and  constraint  (b)  for  choosing 
the  X  's  is  demonstrated  in  example  2.2. 
EXAMPLE  2.2 


Consider  the  circuit  of  example    2  ,  i  , 

1 


(2) 


1 


(3) 


1 
■nrnrtpsy 


(4) 


1/2  £      =S 


^2 


Fig.  2.3 


It  is  desired  to  keep  the  transfer  impedance  as  well  as  the  admittance 
connected  between  nodes  2  and  3  constant. 

From  Fig.  2.2  it  can  be  seen  that  since  x  =  0  corresponds  to 
point  E,  to  stay  in  the  area  of  realizability  we  have  to  choose  X      and 
X       such  that  y??  increases  and  y       decreases  with  x.    Hence  X  ~~ 
must  be  negative  andX~9  must  be  positive,  i.e.  X99  =  -X~~  =  constant 
the  value  of  the  constant  depends  on  the  desired  rate  of  change  of  the 
circuit  elements .    Since  at  the  moment  we  are  not  concerned  with  rate 
of  change  of  elements  we  choose 

V22    =    "y33    =    l 
The  [G1]  ,  [C']and  [  L'      ]  matrices  would  then  become: 
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[G']  = 


-e 


2x 


0 

0 

0 

0 

-2x 

0 

[C']  = 


2x 


-1 


-1 


n      "2X         _ 

3e  0 


[L']  = 


0 

0 

0 

0 

0 

2x 

e 

-1 

0 

0 

.-1 

A      "2X 

4e 

-X 

-e 

0 

0 

-x 

-e 

1 

The  new  circuit  elements  become: 
x 


g'io  "  2"e 

X 

g12   =  e 

2x        x 
9  20    =    e       "e 


9'30    =    e 
g'40    "    l 


-2x 


c' 


c' 


10 


20 


=    2 


2x      . 
e       -  1 


23 


=    1 


c' 


c' 


30 


40 


3e"2x  -  1 


=    2 


20 


23 


L' 


L' 
L' 


30 


34 
40 


2x 
l/(eZX  -  1) 


=    1 

=    l/(4e"2x-e~X  -  1) 

x 
=    e 

=    1/(1  -  e"X) 
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The  variation  of  the  conductances,  capacitances  and  the  inverse  in- 
ductances with  the  independent  variable  x  are  shown  in  Figs.  2.4,  2.5 
and  2.6  respectively. 
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000 


Fig.  2.4     Variation  of  Conductances  with  x,  Example  2.1 
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C30 
Fig.  2.5     Variation  of  Capacitances  with  x,  Example  2.1 
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000 


Fig.  2.6     Variation  of  Inverse  Inductances  with  x,  Example  2.1 
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The  General  Transformation 


3 . 1   Introduction 

In  Chapter  2  a  diagonal  transformation  matrix  was  used  and  it  was 
shown  that  the  realizability  hyper- volume ,  when  it  exists,  can  be 
determined  by  the  manner  described  in  that  chapter.    The  diagonal  trans- 
formation matrix  however  is  limited  in  application.    For  example  it  is 
obvious  from  Equation  (2-8)  that  if  no  conductance  existed  between  nodes 
i  and  j  in  the  original  network,  none  of  the  equivalent  networks  gen- 
erated will  have  a  conductance  between  nodes  i  and  j .    But  the  sim- 
plicity of  the  transformation  (2-3),  linearity  of  the  realizability 
equations  (2-10),  and  the  ease  of  applying  constraints  on  the  transform- 
ed networks  warranted  the  use  of  the  diagonal  transformation. 

In  this  chapter  the  full  n  x  n  transformation  matrix  [T]  is  used. 
The  equations  for  the  elements  of  the  transformed  networks  are  obtained. 
The  realizability  conditions  and  the  allowable  range  for  the  elements  of 
the  transformation  matrix  [T]   so  as  to  yield  to  equivalent  networks 
with  positive  R,  L,  C  circuit  elements  is  then  discussed. 
3  .2  The  General  Transformation  and  Realizability  Conditions 

The  transformation  matrix    T  is  chosen  to  be 


W- 


11 


21 


y 


y 


12 


22 


In 


2n 


*nl 


nn 


(3-1) 


A  typical  element  of  the  conductance,  capacitance,  and  inverse 
inductance  matrices  [G1]  ,  [C]  and  [  L1      ]  of  the  equivalent  networks 
produced  will  then  be  of  the  form: 
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g'ij  ■  ,\  ,  wilc'lk  (3-2) 

£  ,k=l 

Since  the  original  network  is  made  up  of  bilateral  elements,  the 
[G  ]  ,   [C  ]  and  [L      ]  matrices  are  symmetric.    Hence  Equation  (3-2) 
can  be  written  as: 

n  n 

g\.    =     2     g, .  y.,  y  ..   +     E         g       [y.    y.     +y.    y.     ]  (3-3) 

y  U  ykk'ik     jk  pq     'ip'jq         lq     ]P 

q=p+l 

Where  i  =  j,  Equations  (3-3)  yield  the  diagonal  elements  of  the  matrices 

[G'  ]  ,   [C   ]  and  [L,_1]  .    Thus: 

n  n 

g'..    =     £     g,  .  y.,      +2      L        g      y .    y .  (3-4) 

y  ii         k=1    ykk 'ik  ypq'ip'iq 

q=p+l 

We  recall  from  Chapter  1  that  realizability  with  passive  components 

required  that  the  off  diagonal  elements  of  the  [G*  3  ,   [  C  ]  and  [  L'      ] 

be  zero  or  negative  and  that  the  sum  of  terms  in  any  row  of  the  above 

matrices  be  zero  or  positive.    Thus  realizability  with  passive  elements 

in  the  case  of  the  full  transformation  matrix  [  I"  ]  requires  that 

g'..       ^     0  for  i  ^   j  (3rS) 

n 
g'.      =    I    g'       ^     o  (3-6) 

10    j-1    1J 

For  an  n  port  network  each  of  the  matrices  [  G'  ]  ,  [  C'  ]  and 
[  L'  ]  is  a  symmetrical  n  x  n  matrix.  Hence  we  have  in  general 
3n(n-l)/2  inequalities  of  the  form  (3-5)  and  3n  inequalities  of  the  form 


(3-6)  to  be  satisfied.     Since  the  transformation  matrix  [T]   need  not  be 

2 
symmetrical  it  contains  n     elements.    The  invariance  of  the  input  or 

transfer  impedance  requires  that  one  or  two  columns  of  the  transformation 
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matrix  [  T  ]  be  unit  vectors  (i.e.  zeros  everywhere  except  on  the  diag- 
onal position).    Hence  if  the  invariance  of  the  transfer  impedance  is 

required  a  total  of  3n(n+l)/2  inequalities  have  to  be  satisfied  by 

2 
(n  -2n)  =  n(n-2)  y  's  if  the  transformed  network  is  to  be  realizable  with 

passive  components. 

Examination  of  (3-3),   (3-4),   (3-5)  and  (3-6)  shows  that  the  in- 
equalities (3-5)  and  (3-6)  which  have  to  be  satisfied  are  non-linear 
functions  of  the  y  's.    Hence  the  method  used  in  Chapter  2  for  the  diag- 
onal transformation  matrix  is  not  applicable  here,  since  the  equations  to 
be  solved  are  not  simultaneous  linear  equations.    To  be  able  to  apply 
the  method  of  Chapter  2  to  the  transformation  of  (3-1)  used  in  this  chap- 
ter, we  will  linearize  the  inequalities  of  (3-5)  and  (3-6)  as  shown  in 
the  next  section.    It  must  be  noted  that  in  the  case  of  transformation 
with  the  full  transformation  matrix  of  (3-1), in  contrast  to  the  case  of 
diagonal  transformation  matrix ,the  y  's  can  take  both  positive  and  neg- 
ative values . 
3  .3  Linearization 

We  know  from  the  discussion  in  Chapter  2  that  the  inequalities 
(3-5)  and  (3-6)  are  satisfied  when  y. .    =    1    and   y,.    =    0  f or  i  ^  j . 

•xThis  will  result  in  the  transformation  matrix  [T]  to  be  the 
identity  matrix  [  I  ]  and  the  transformation  will  yield  the  original  net- 
work. 

We  can  linearize  the  inequalities  (3-5)  and  (3-6)  by  expanding  the 
matrix  [Y1]  in  a  Taylor  series  expansion  about  the  pointy.,    =    1  and 
y..    =    0       i    /j    as  shown  below: 

From  (1-3)  we  have 

[Y'l 
or  from  (l-6a)  we  have 


[Y']   =    [T][Y  ][r]T  (3-7) 


LG']   =    [r][G  ][r]T  (3-8) 
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Hence 


[Ac]  =  [Ar][G  ][r]T  + [r][G  ][Ar]T 


:.[AG']  =  /     £ 


i/j=l  ij 


ay.,  ij 


n 
x        I 


— r A  y ., 


(3-9) 


But 


o[r] 
by  . 


u 


Ayy   - 


[ 


Ay., 
ij 


Hence 


n 
L 

i=l 
3=1 


~S Ay  .. 


Ay 


Ay 


11 


21 


Ay  , 

'  nl 


Ay 


Ay 


12 


22 


Ay 


In 


Ay 


nn 


=LAy] 


(3-10) 


Since  we  are  expanding  [G']  with  a  Taylor  series  about  the  point 


y . .  =  1   #  y . .  =  0  i^i       we  can  write 

li  ij 


v..    =    1+6.. 
n  n 

v..    =    0+6.. 


(3-11) 
(3-12) 
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where  6..  and  6..  are  incremental  change  ony,,  andy,.   .i.e. 
11  ij  'ii  ij 

6..    =    Ay.,  and    6..    =   Ay  .    and  the  transformation  of  (3-1)  with 
n  n  ij  ij 

6..    =    6..    =    0    yields  the  original  network .    The  origin  of  the  6-space 
corresponds  to  [  T  ]   being  the  identity  matrix. 

Substituting  (3-11)  and  (3-12)  in  3-9  we  can  write 

[6]  =    [Ay]  =    I  \[T  ]       Ay..  (3-13) 

U-l         ^ij  1J 

Substituting  (3-13)  in  (3-9)  we  obtain 

[AG']  =    [6][G  ][T]T+  [T][G  ][6  ]T  (3-14) 

Expanding  [G1]  with  the  Taylor  series  around  [G  ]  yields 

[G']  -    [G']  |       +    [AG']  |  (3-15) 

o  o 

where  [G1]  I      and  [AG1]  I      denotes  the  matrices  of  [G']and 
o  o 

[  AG '  ]  where    y..    =    1    andy..    =    0        i^j    or  where  6  . .  =  6 . .  =  0 . 
ii  ij  ii         ij 


From  (3-8) 


EG']  |Q   =    [G  ]  (3-16) 


Since  [T]  I      =[r]T        =[U]  where    [U]  denotes  the  identity 
o  o 

matrix  then  from  (3-14)  we  have 


[AG1]  |q    =    [6][G  ]   +    [G  ][6]T  (3-17) 

Substituting  (3-17)  and  (3-16)  in  (3-15)  yields 


[G']=    [G  3  +    [6][G]+    [G][6]T  (3-18) 

Hence  a  typical  term  of  the  [G1]  ,  [C]  and  [  L1      ]  matrices  would 


be 


n 


""u-'ii*,^   (gki6ik  +9ikV  (3-19) 

where  only  the  first  two  terms  of  the  Taylor  series  are  taken. 
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Equation  (3-19)  will  enable  us  to  linearize  the  inequalities  of 
(3-5)  and  (3-6)  and  hence  apply  the  method  of  Chapter  2  to  obtain  the 
area  of  realizability  and  the  range  of  y  's  for  which  the  transformation 
of  (3-1)  yields  equivalent  networks  with  passive  elements.    In  essence 
this  procedure  identifies  a  small  region  in  6-  space,  any  point  of  which 
corresponds  to  a  set  of  y  's  for  which  the  transformation  of(3-l)may  be 
applied  to  yield  equivalent  networks  with  passive  elements.    Once 
this  region  is  defined  a  searching  function  such  as  the  exponential 
searching  function  of  Chapter  2  can  be  used  to  generate  equivalent 
networks  whose  circuit  elements  are  functions  of  an  independent  variable 
such  as  x  and  the  complete  inequalities  of  (3-5)  and  (3-6)  investigated 
to  yield  equivalent  networks  with  passive  elements. 
3  .4  The  Exponential  Searching  Function 

We  use  an  exponential  searching  function  for  the  elements  of  the 
transformation  matrix  [  I  ]  in  the  following  form: 


X..  x 


y..    =   e  "  (3-20) 

X..  x 
Vy    -    1  "  e  1J  i;«j  (3-21) 

Approximating  (3-2  0)  and  (3-21)  by  the  first  two  terms  of  their  Taylor 
series  expansion  and  comparing  with  Equations  (3-11)  and  (3-12)  respec- 
tively we  obtain 

X..  x    =    6..  (3-22) 

n  n 

-X..  x    =    6..  (3-23) 

U  ij 

Equations  (3-22)  and  (3-23)  in  conjunction  with  the  linearized  form 
of  Equations  (3-5)  and  (3-6)  can  be  used  to  choose  values  of  the  X  's 
and  the  range  of  x,  as  best  illustrated  by  the  following  example. 
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EXAMPLE  3.1 

Consider  the  RC  network  of  Fig.  3.1  . 


(1) 


-n/V\AAAA/v 


— c 


Fig.  3.1 


(2) 


It  is  desired  to  find  the  equivalent  networks  containing  passive  elements 
which  have  the  same  input  impedance  as  the  above  network. 
The  [G  ]  and  [C  ]  matrices  of  the  above  network  are: 


[C  ]  = 


c 


0 


1 


0 

c, 


[G  ]=g 


1 
-1 


-1 
1 


Starting  with  the  diagonal  transformation  matrix 

1  0 


[T]  = 


0 


y 


22 


The  [G1]   and  [C]  matrices  will  be 
1 


[G']=g 


22 


"y22  V22 


CO 


C2V22 


The  realizability  conditions  are 


1  -y 


22 


^      0 


"1+y22     &      ° 
C2y222        ^      0 
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It  is  obvious  that  the  only  value  oiy      which  satisfies  the  above 
equations  is  y       =  1  which  means  that  the  transformation  matrix  is  the 
identity  matrix.    Hence  the  diagonal  form  of  the  transformation  matrix 
will  not  generate  any  equivalent  network. 

Now  we  use  the  full  transformation  matrix 


[r]  = 


1 


12 


22 


The  [G1]  and  [C]  matrices  will  be 


[G']=g 


l-2y12+y12 


y22  V12  "  V22 


y22yl2      y22 


22 


[C']  = 


r 


C1+C2*E2 


C2V12  y22 


C2V12y22 


2 
C2y22 


The  realizability  conditions  (3-5)  and  (3-6)  become: 
(a)    g'l0  =  (1-2v12+yi22+yl2>'22->/22)g 


{b)    g,20==(Y22yl2-Y22+>/22)g      * 
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(c)  s'12  =  (>'22y12  -y22)g    «      P 


(d)    C'lO  =  Cl+C2(y12y22+yl22)      *        ° 


(e)    C'20  =  C2(y12y22  +y22   >       *        ° 


(f)     C'l2  =  C2^12y22      < 


The  linearized  form  (3-19)  of  the  above  equations  are 
(a)    "  612  "622       *     ° 


(b)    612+622       »      0 


(c)    l+622-612       »      0 


(d)    C1+C2612       i      0 


(e)    C2(612  +  2  622  +  l)       £      0 


(f)   c2e12     <    0 


Each  of  the  above  inequalities  when  satisfied  with  equal  signs  will  be 
presented  by  a  line  on6,„  ,  6    .  space  as  shown  in  Fig.  3.2.    The 
linearized  realizability  region  reduces  to  the  portion  of  the  line 
699  +  6       =  0    in  the  second  quadrant  of  the  6  space  shown  on  Fig.  3.2 

Using  the  exponential  search  function,  Equations  (3-20)and (3-21) 
become 


X  22  x 


v22    =    e 


,         X12X 
^12    =    l~e 
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Equations  (3-22)  and  (3-23)become 


X22X    =    622 


"X12X    =    612 


From  the  conditions  of  realizability  for  this  case  we  have  6~9  =  -6   9  = 

constant.    The  value  of  the  constant  determines  the  rate  of  change  of 

the  numerical  values  of  the  circuit  element.    Since  we  are  not  concerned 

at  this  moment  about  the  rate  of  change  of  elements  we  arbitrarily 

choose  X92  =  X  ?  =  1    for  x  >  0 . 

Hence 

x 
y22=e 

1  X 

yi2=1-e 
Substituting  the  above  values  in  the  expressions  for  g's  and  c's  yields 
9'10-0  C'10  =  C1  +  C2(l-eX) 

^20  =  °  C'20  =  C2eX 

g'12  =  -ge2x  C'12  =  C2ex(l-eX) 


From  the  above  equations  it  can  be  seen  that  the  range  of  x  for  which  the 
realizability  conditions  are  satisfied  is  given  by  the  two  equations 

Cl  +C2(1  -  eX)       £      0 
C2eX(l  -  eX)  <      0 


Hence 


ci 


1       ^      e         ^      1  + 


C2 


c    +  c„ 

1  2 

or  0       <      x      -$       I  n    ■ 


C2 
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The  equivalent  network  is  shown  in  Fig.  3.3  below  and  the  variation  of 
the  network  elements  with  x  are  shown  in  Fig.  3.4  for  the  case  of  g  =  1 , 

1=C2 


c-    "C„"l.  C„eX(eX-l) 


■II- 


/^s^/s/NAA/^A- 


ge 


2x 


x, 


Cj+C^l-e  ) 


-rC2e 


Fig.  3.3 
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Fig.  3.4     Variation  of  Circuit  Elements  of  Example  3.1 
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Conclusion  and  Suggestions  for  Further  Research 
4 . 1    Conclusion 

In  this  paper  the  congruent  transformation  is  applied  to  the  admit- 
tance matrix  of  an  n+1  node  common  datum  network  to  generate  equivalent 
networks  which  have  the  same  input  or  transfer  admittance.    It  is  shown 

that  a  transformation  matrix  whose  elements  are  exponential  functions 

Xx 

of  a  real  variable  x,  such  as  e      ,  can  be  used  to  generate  a  group,  of 

equivalent  networks  whose  circuit  elements  are  functions  of  the  variable 
x.    The  range  of  variation  of  the  elements  of  the  transformation  matrix 
required  to  generate  equivalent  networks  with  passive  circuit  elements 
is  found  by  solving  a  number  of  sets  of  equations  called  the  realizability 
equations.  If  a  diagonal  transformation  matrix  is  used  these  equations 
are  linear  but  when  a  full  transformation  matrix  is  used  these  equations 
are  non-linear,  and  have  to  be  linearized.    If  a  diagonal  transformation 
matrix  whose  elements  are  exponential  functions  of  the  real  variable  x 
is  used,  then  certain  constraints  can  be  placed  upon  the  generated 

equivalent  networks  by  choosing  certain  relations  between  the  constant 

Xx 
factors  (  X  's)  in  the  exponent  of  the  exponential  function  e 

4  .2   Suggestions  for  Further  Research 

It  is  believed  that  the  subject  of  equivalent  networks  has  a  great 

potential  for  application  in  the  field  of  network  synthesis  and  analysis. 

The  following  are  some  ideas  for  further  research  on  this  subject. 

(a)  It  may  be  possible  to  find  the  general  form  of  a  transformation 
which  would  take  a  given  form  of  a  synthesised  network  to  another  form., 
For  example  it  may  be  possible  to  find  a  transformation  matrix  [I"*]  which 
transforms  the  Foster  form  of  a  synthesised  network  to  Cauer's  form  for 
the  same  input  impedance. 

(b)  The  effect  of  the  presence  of  mutual  inductance  in  a  network 
on  the  transformation  can  be  investigated  with  the  aim  of  transforming 

a  network  with  mutual  inductance  to  a  network  without  mutual  inductance. 

(c)  Finally,  from  curves  of  Fig.  2.4,  2.5  and  2  . 6  it  can  be  seen 
that  if  the  values  of  the  components  of  the  network  vary  according  to 
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the  curves  shown,  then  the  transfer  impedance  and  input  impedance 
remain  invariant,  i.e.  the  multiparameter  sensitivity  of  the  input  or 
transfer  impedance  is  then  zero.    Hence,  multiparameter  sensitivity 
can  be  studied  from  the  point  of  view  of  equivalent  networks. 
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Appendix  I 
Invariance  of  Input  Impedance 
It  was  stated  in  Chapter  1  that  if  the  column  i  of  the  transformation 
matrix  [  l  ]  is  a  unit  vector  (zero  everywhere  except  for  a  one  at  the  1th 
row)  then  the  input  impedance  at  port  i,  that  is  Z. ,  remains  invariant 
under  the  transformation.    In  this  appendix  a  proof  of  this  fact  based  on 
matrix  portioning  is  presented. 

Consider  an  n+1  node  network  with  a  current  source  connected  be- 
tween each  of  the  n  nodes  and  the  datum  node  as  shown  in  Fig.  1-1 . 


Fig.  1-1 

The  equilibrium  equation  of  the  network  can  be  written  in  matrix 

form  as 

[I]-[Y][V-]  (|-l) 

with  [  I  ]  ,  [  Y  ]  and  [  V  ]  as  defined  in  Chapter  1 . 

Expanding  the  matrix  Equation  (1-1)  and  using  Crammer  rule,  the 
voltage  V    between  node  i  and  datum  is  given  by 

where: 

A     =     the  determinent  of  the  admittance  matrix  [  Y  ] 
A. .  =     the  minor  formed  from  Aby  deleting  row  i  and  column  j 
If  the  network  is  excited  only  at  node  i  by  a  current  source  I.  then; 
I .  =  0         for  i  =  1    .    .    .    .  n         i  ?  j 


i 
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Hence 
and 


A.. 

V.  -  -IT-  I. 
1  A        1 


(1-3) 


A. 


Z.  = 

i 


li 


(1-4) 


The  invariance  of  the  input  impedance  at  port  i  under  the  congruent 
transformation 


[y]= [r][Y ][r] 


(1-5) 


with  column  i  of  [  I  ]  being  a  unit  vector  would  be  proved  if  we  can 
show  that  the  ratio    A..  /A  of  Equation  (1-4)  remains  invariant  under 
transformation  (1-5). 

For  simplicity  of  proof  we  consider  the  invariance  of  the  input  im- 
pedance at  node  1.    Since  by  renumbering  the  nodes  we  can  call  any  node 
node  1 ,  there  is  no  loss  of  generality  and  the  proof  would  be  acceptable 
for  the  invariance  of  the  input  impedance  at  any  port. 

Consider  a  real  non-singular  transformation  matrix  [T]  where 
column  1  is  a  unit  vector 


[r> 


i 


12 


22 


In 


2n 


0         y 


n2 


.     y 


nn 


(1-6) 


The  admittance  matrix  [  Y' ]   of  the  transformed  network  is  given  by 


[y-]=  [iHYHr]1 

where  [  Y  ]  and  [Y1]  are  as  defined  in  Chapter  1 . 


(1-7) 
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Since  the  determinant  of  the  product  of  matrices  is  equal  to  the 

product  of  the  determinants  [!"']  then: 


A1  =   |r 


(i-8) 


where    A'  =  determinant  of  [  Y' ] 
and   |  r  |    =  determinant  of  [  I  ] 

We  can  partltfcn  the  [  T  ]  ,  [  T  ]      ,  [  Y  ]  ,  [  Y'  ]  matrices  as  shown 
below.* 


11 


y' 


In 


X' 


y-  I 

nl    ! 


11 
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1      1    V 
I    y12 

_J 

°  ; 


o     i 


In 


0 

•                     • 

• 

0 

12  1 

T 
a 

In 

Performing  the  matrix  operation  we  obtain 

[X']=  Co][X  ][a]T 
Hence  the  admittance  matrix  of  the  new  network,  is  given  by: 


(1-9) 


(I- 10) 


y. .    i  ^  j    =    the  admittance  connected  between  nodes  i  and  j 

y..    =    the  sum  of  all  admittances  connected  to  node  i 
li 
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[Y']  = 


11 


In 


i    [«][X  ][<*]T 


y' 


nl 


(1-11) 


The  input  impedance  at  node  i  of  the  new  network  can  be  written 
from  Equation  (1-4)  to  be 


A* 


Z'l  = 


11 


det.   (  [gj  [X  ][g]x) 
A' 


(1-12) 


Since  determinant  of  product  of  matrices  is  equal  to  the  product  of  the 
determinants  [  7  ]  we  can  write 


det  (  [a][X][a]T)  =  |a|2|X  1=  |a|2  A 


11 


(1-13) 


where    |  a  |   =    determinant  of  [  a  ]  . 

Since  the  1st  column  of  [  I  ]  is  a  unit  vector 

|r|  =    lor  I 

Substituting  (1-14)  in  (1-13)  we  obtain 


i 


det  (  Lo][a][a]T)  =  If  |2  A 


11 


Substituting  (1-8)  and  (1-15)  in  (1-12)  yields 


z'i  = 

UI2An 

|r|2  a 

Hence " 

■ 

*1 


z"l  =  zl 


(1-14) 


(1-15) 


(1-16) 
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Appendix  II 
Computer  Program  for  Example  2.1 


..JOB0184F.ARDALAN 

PROGRAM  EQTRANF1 

DIMENSIONX1500) .010 (500) .012 (500) ,G20< 500) .G30( 500) »G40<500) ,C10( 
1500 ) »C20( 500 )»C23( 500 )»C30( 500) .C40C500) »P20(500) »P23<500) »P30( 50 
2)»P34(500)»P40(500).ITITLE(12) .GA2 ( 500 ) »GA3( 500 ) 

GM2=1. 

GM3«-1. 

X(1)*0. 

A»0.8 

DO  10  1-1,500 

G12(I)»EXPF(GM2*X< I) ) 

G10(I)«2.-G12(I) 

G20( I)=EXPF(2.*GM2*X( I ) )-G12( I ) 

G30<I)=EXPF(2«*GM3*X(  I  )  ) 

G40( I)  =  l. 

C10< 1)»2. 

C23(I)=EXPF( (GM2+GM3)*X( I )  ) 

C2QU)=EXPF(2«*GM2*X<  I))-C23(I  ) 

C30(I)»3**EXPF(2«*GM3*X(I> >-C23(I> 

C40U)»2« 

P23(I)»EXPF( (GM2+GM3)*Xll) ) 

P20(I )=EXPF(2.*GM2*X< I))-P23(I) 

P34< I)=EXPF(GM3*X( I) ) 

P30<I)»4.*EXPF(2.*GM3*Xm  J-P23  (  I  )-P34<  I) 

P40(I )*1.-P34( I ) 

GA2( 1)«EXPF(GM2*X< I) ) 

GA3(I)»EXPF(GM3*X< I) ) 

X(I+1)=X( I)+(A/500«) 
10  CONTINUE 

DO  22    1*1*12 

22  ITITLE(I)»12 
ITITLE(1)*8HARDALAN 
LABLE=4HG10 

CALL  DRAW ( 5G0.X. G10. l.Q.LABLE. IT ITLE.0. 0.0.0.0.0.5. 5.0. LAST) 
DO  23  1=1.12 

23  ITITLE(I)=12 
ITITLE( 1)«8HARDALAN 
LABLE=4HG12 

CALL  DRAW(500»X»G12»2»0»LA8LE.lTlTLE»O»O»O»0t0»0,5*5.0»LAST) 
DO  24  1=1.12 

24  ITITLE( I >  =12 
ITITLE(1)»8HARDALAN 
LABLE«4HG20 

CALL  DRAW(50O»X.G20.2»0.LABLE. I T I TLE .0.0 »0 .0 .0 .0. 5 . 5 »0 . LAST ) 
DO  25  I =1.12 

25  ITITLE(I)=12 
ITITLE( 1)=8HARDALAN 
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LABLE=4HG30 

CALL  DRAW(500»X.G30.2»0.LABLE.lTITLE.0.0.0.0,0»0»5,5»0,LAST) 

DO  26  1=1,12 

26  ITITLE(I)=12 
ITITLE(1)"SHARDALAN 
LABLE=4HG4Q 

CALL  DR AW ( 5 00 » X • GAO , 3 » 0 » LABLE , I T I TLE .0 . 0 .0 .0 .0 • 0 ♦ 5 • 5 .0  » LAST ) 
DO  27  1=1*12 

27  ITITLE(I)=12 
ITITLE( 1>»8HARDALAN 
LABLE=4HC10 

CALL  DRAW<500.X»C10.1»0.LABLE»ITITLE.0,0.0.Q.0.0.5.5»0»LAST) 
DO  28  I=l»12 

28  ITITLE(I)=12 
ITITLE(1)»8HARDALAN 
LA8LE=4HC23 

CALL  DRAW < 500 .X.C23. 2 .O.LABLE. IT ITLE .0.0 .0.0.0 .0.5. 5.0 .LAST) 
DO  29  1=1*12 

29  ITITLE(I>=12 
ITITLE(1)=8HARDALAN 
LABLE=4HC20 

CALL  DRAW < 500 .X.C20.2 .O.LABLE .ITITLE .0.0. 0.0.0.0. 5.5.0 .LAST) 
DO  30  1=1.12 

30  ITITLE(I)=12 
ITITLE(1)=8HARDALAN 
LABLE=4HC30 

CALL  DRAW(5Q0»X.C30,2.0»LABLE.lTITLE»0.0»0.0.0.0.5»5.0»LAST) 
DO  31  1=1.12 

31  ITITLE(I)=12 
ITITLE( 1)=8HARDALAN 
LABLE=4HC40 

CALL  DRAW ( 500 .X.C40. 3 .O.LABLE. I T I TLE .0.0.0 »0.0 .0, 5. 5 .O.LAST ) 
DO  32  1=1.12 

32  ITITLE(I)«12 
ITITLE(1)=8HARDALAN 
LABLE=4HL23 

CALL  DRAW( 500 .X.P23.1. O.LABLE. U I TLE .0.0. 0.0 .0.0. 5. 5 .O.LAST ) 
DO  33  1=1.12 

33  ITITLE(I)=12 
ITITLE(1)»8HARDALAN 
LABLE=4HL20 

CALL  DRAW( 500. X »P20. 2 »0 »LABLE» IT ITLE» 0.0.0.0, 0.0.5. 5 »0» LAST) 
DO  3h    1=1.12 

34  ITITLE(I)=12 
ITITLE(1)=8HARDALAN 
LABLE=4HL34 

CALL  DRAW( 500. X.P34.2.0. LABLE, ITITLE.O. 0.0. 0.0. 0.5.5.0. LAST) 
DO  35  1=1.12 

35  ITITLE(I)=12 
ITITLE(1)=8HARDALAN 
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LABLE=4HL30 

CALL  DRAW<500tX,P30»2»0»LABLE.lTITLE»0,0»0»0.0.0»5t5»0»LAST) 
DO  36  I«l#12 
36  ITITLE(I)=12 

IT1TLE( 1)*8HARDALAN 

LABLE=4HL40 

CALL  DRAW(5  00.X,P<*0»3.0.LA8LE, I  T  I  T LE »0, 0 tO »0 ,0 .0, 5 • 5 »0 ♦ LAST ) 

END 

END 
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Appendix  III 
Computer  Program  for  Example  3.1 

..JOB0184F.ARDALAN 

PROGRAM   NONLINER 

DIMENSION   X(500).G1Q<500) .612(500) .G20< 500 ) ,C10( 500 ) »C12 ( 500 ) .C20< 
1500).ITITLE(12) 

A=l. 

GM12=1. 

GM22«1. 

G»l« 

Cl*l. 

C2  =  l. 

X(1)«0. 

00    10    1=1,500 

G12(l)=G*EXPF< (GM22+GM12)*X(I)) 

C12( I)=C2*(EXPF(GM22*X(I) ) *( EXPF(GM12*X ( 1) )-l.) ) 

C10(I)*C1+C2*(EXPF(GM22*XU) )*( 1,-EXPF ( GM12*X< I) ) )+ ( 1 .-EXPF ( GM12*X 
1(1))  )**2) 

C20(I)=-C12( I )+C2*EXPF(2.*GM22*X( I ) ) 

PRINT  15»X( I ) »G12( I) »C10( I)»C1Z( I) »C20(  I) 
15  FORMAT ( 1X.5F10.6) 

X( I+1)=X( I)+ (A/500.) 
10  CONTINUE 

DO  23  1-1.12 

23  ITITLE(I)=12 
ITITLE(1)*8HARDALAN 
LABLE=4HG12 

CALL  DRAW ( 500.X.G12. 1.0.LABLE, IT 1TLE.0.Q .0.0 .0.0.5 .5.0 .LAST  ) 
DO  24  1=1.12 

24  ITITLE(1)»12 
ITITLE(1)=8HARDALAN 
LABLE=4HC10 

CALL  DRAW( 500 .X.C10. 2 .O.LABLE. IT ITLE .0.0.0. 0.0.0. 5.5  .O.LAST) 
DO  25  1*1.12 

25  ITITLEC I )*12 
ITITLE(1)=8HARDALAN 
LABLE=4HC12 

CALL  DRAW(  500  .X.C12.-2  .O.LABLE.  I  TITLE. 0.0. 0.0  .0.0.  5.  5. O.LAST) 
DO  26  1=1.12 

26  ITITLE(I)=12 
ITITLE(1)=8HARDALAN 
LABLE=4HC20 

CALL  DRAW(500.X.C20»3.0.LABLE.ITITLE.0.0.0.0.0»0.5.5»0.LAST) 

END 

END 
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